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Abstract 

We discuss conditions under which expectation values computed 
from a complex Langevin process Z will converge to integral averages 
over a given complex valued weight function. The difficulties in prov- 
ing a general result are pointed out. For complex valued polynomial 
actions, it is shown that for a process converging to a strongly station- 
ary process one gets the correct answer for averages of polynomials if 
Cr(fc) = E{e^^^^'^^) satisfies certain conditions. If these conditions are 
not satisfied, then the stochastic process is not necessarily described 
by a complex Fokker Planck equation. The result is illustrated with 
the exactly solvable complex frequency harmonic oscillator. 



* Supported by Fonds zur Forderung der Wissenschaftlichen Forschung in Osterreich, 
project P7849. 
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1 Introduction 



It is well known that a straightforward application of standard simulation 
techniques like the Metropolis et al. algorithm |l| will fail when they are ap- 
plied to problems with a complex action or Hamiltonian. This is due to the 
fact that there is no direct probabilistic interpretation of a distribution func- 
tion of the form when the action S is complex valued. As an alternative 
to such algorithms, the complex Langevin (CL) method was first proposed 
by Klauder , and subsequently studied by many authors. The central idea 
of CL is based on the fact that for a Langevin equation there is no formal re- 
striction to a real valued drift term. The use of a complex drift term provides 
CL with a genuine advantage over other methods; namely, since CL uses the 
entire complex action S to define a stochastic process, it can in principle 
converge directly to the desired distribution. This potential for circumvent- 
ing the well known "sign problems" of other standard algorithms is just one 
reason CL continues to be a subject of great interest. 

Quite independent of its utihty as a numerical technique is the interesting 
fact that under certain conditions a system governed by a complex Hamil- 
tonian can nevertheless be given a probabilistic interpretation. Indeed, the 
name complex Langevin may be slightly misleading, since the Langevin equa- 
tions really describe a real diffusion process in twice as many dimensions. 

Unfortunately, there is currently no complete theory of the CL method. 
Many conditions under which a real Langevin process can be shown to con- 
verge to a given distribution are not satisfied for a general CL process. Fur- 
thermore, from the point of view of numerical simulations, in some cases CL 
is known to converge to the wrong results For simple actions, this 

truant behavior can be corrected by an appropriate choice of kernel in the 
Langevin equation 0, but for more general systems, in particular lattice 
fermion models of current interest, it is far from clear which choice of kernel 
is required. 

The purpose of this paper is to explore in a rigorous fashion the condi- 
tions under which the CL process correctly simulates a given system with a 
complex Hamiltonian defined on a Euclidean space a; G M". While a general 
theorem is still lacking, we will demonstrate a set of sufficient conditions for 
polynomial actions. 
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2 The Process 



Throughout the paper we will assume that the system is described by vari- 
ables X G with some complex action or Hamiltonian S : MP' C, where 
KeS is bounded from below. For simplicity the discussion will be restricted to 
one degree of freedom, since the following statements allow for an immediate 
generalization to M". The quantities of physical interest are of the form 

(gi^)) = ^ J^9{x)e-'^-Ux, (2.1) 
f e-sWrfx , (2.2) 



assuming that for the partition function Af we have < lA/"] < oo. Both 
g{z) and S{z) are assumed to be analytic in C. Since S{z) is analytic, this 



provides a local Lipschitz condition for the Langevin equation p.8| ; that is, ^ 
has a unique local solution that is defined up to a random explosion time 
For real actions G M we define a process {X{t),t > 0} by the Langevin 
equation 

dX{T) = F{X{T))dT + dW{T) , (2.3) 

with the drift term 

, IdS(x) , 

where Wir) is a standard Wiener process with zero mean and covariance 

< W(T^)W{r2) >= min(ri, r^). (2.5) 

The probability density f{x, r) for such a process will converge pointwise to 
the desired distribution function. That is 

lim f{x,T) = f{x) a.e. , (2.6) 

with 

f{x) = le-^(^) . (2.7) 

As already mentioned in the introduction, for the complex case one can 
formally construct a Langevin equation 

dZ{T) = F{Z{T))dT + dW{T) , (2.8) 
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with the drift term 

.(-)^-i^. (2.9, 

As above, W{t) is a standard Wiener process. This is actually a two dimen- 
sional process of the form 

dX{T) = G{X{t), Y{T))dT + dW{T) (2.10) 

dY{T) = H{X{t), Y{T))dT (2.11) 
with S{z) = u{x, y) + iw(x, y) 

G(.,,).-i^, H(.,y)^\^. (2.12) 

Note that the equation for dY has a zero diffusion coefficient, but is never- 
theless a stochastic equation through its dependence on X . 



The process {(X(r), F(r)), r > 0} as defined by equation has a dis- 
tribution density f{x,y,T). There are now two crucial questions concerning 
this process. The first question is whether this so defined process converges 
in distribution at all to some {X,Y), 

lirn^ f {x, y, t) = f{x,y), a.e. . (2.13) 
The second question has to do with the problem of whether f{x, y) satisfies 
E{g{X + zY)) = J^^ g{x + iy)f{x, y)dxdy = ^ j^g{x)e-'^-Mx. (2.14) 

This equation contains the essence of complex Langevin. It tells us that if 
the process as defined above has converged in distribution, we might be able 
to calculate {g{x)) by an equivalent probabilistic system of twice as many 
dimensions. 

Before we investigate equation |2.14| it is necessary to discuss the asymp- 



totic behavior of X(t),Y(t) for r oo. To do this we will examine the 
equivalent Fokker-Planck (F-P) equation. First note that since the diffu- 
sion and drift coefficients are independent of r, the process (X(r), Y{t)) is a 
homogeneous diffusion process. For twice continuously differentiable distri- 
bution densities with respect to x, y and once with respect to r, there exists 
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a F-P equation . Then the equivalent F-P equation for the above Langevin 
equation is 

dfix,y,T) 



dr 

with 



T/(x,i/,r) (2.15) 



T = -G^{x,y) -G{x,y)—- Hy{x,y) - H{x,y)— + --^ . (2.16) 

As may be seen from the above equation, this case only requires a continuous 
first order derivative with respect to y. Let us first assume that T has a 
unique stationary solution f{x,y) 

Tfix,y) = , (2.17) 

with 

/(a;,2/) >0, a.e. (2.18) 

and 

f{x,y)dxdy = 1 . (2.19) 



One can then use f{x,y) to define an invariant measure 

//(0)= / f{x,y)dxdy (2.20) 



with respect to Lr = exp{QT), r > , which is defined as an operator 
family in L^(M^, dfi) {p = 1,2). The operator Q is obtained by transforming 
T to Li(M^rf/i) by 

Q = e^Te-^ . (2.21) 

Then there is a theorem which tells us that {L^-^r > 0} is a contraction 
semigroup 0, from which it follows that for any function G U'{U?'^dii) , 
p = 1, 2 one has convergence in the strong sense 

s- lim = cs , (2.22) 

where is a constant. Since /i(M^) < oo, any distribution density f{x,y,T) 
converges pointwise to the stationary solution 

lim f{x,y,T) = f{x,y) a.e. . (2.23) 
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On the other hand, if the zero eigenvalue of T is M-fold degenerate then 
there are M ergodic classes. In this case we have 



M 

Jim /(x, r) = ^ Cifi{x, y) a.e. . (2.24) 

Unfortunately, it is still an open question as to what conditions guarantee 
the existence of a stationary solution at all. Note that the CL process defined 
by the real two dimensional equations p.lQ and |2.11| has a singular diffusion 



matrix. Therefore, a general statement on the existence of stationary solu- 
tions can not be made based on the nature of the drift terms. 

Further, T does not fall into the class of hypoelliptic operators of constant 
strength for general S{x + iy), {x,y) E R^. In this general statement 

on the regularity of the solutions can not be made [§], thus one can not 
exclude the possibility that the solutions Tf = exist only in the sense 
of distributions (weak solutions). If this is the case then it might be quite 
difficult to find (construct) the stationary density. 

It is easy to see that a situation like this can occur with the very simple 
example S{x) = cx^, c G , which is now supposed to be solved by 
complex Langevin. The Langevin equation reads: 

dX{T) = -cX{T)dT + dW{T), (2.25) 

dY{T) = -cY{T)dT (2.26) 

The stationary density f{x, y) is then a weak solution and can be formally 
given as 

/(a;,2/)~e-^^'5(y). (2.27) 



3 Polynomial Actions 



Returning to the question of equation ^.14 , we now demonstrate the following 



general result: let S{x) and g{x) be polynomials of degree and M (M < 
— 1) and S{x) such that G iS(M). iS(M) is the Schwartz space of 
functions of rapid decrease. Assume Z{t) converges in distribution to 
a strongly stationary process. Then equation |2.14| holds if there exists a Tq 
such that for all r > Tq 
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1. 



2. The Fourier transform 

h{x, T 

of the expectation value 



< oo for all < n < - 1, /c e 



1 

2^ 



Cr{k)e-''''dk 



(3.1) 



(3.2) 



(3.3) 



satisfies h{x,T) G C^(]R) with respect to x, and h{x,t) G C"^(M) with 
respect to r. Furthermore, x'^~^h{x,T) G -/^-"^(M, rfx) . 



3. 



hm Cr{k) G 5(M) 



(3.4) 



Before proving this result, we note that condition 3 is also a necessary 
condition for equation |2.14| , since, by assumption, G S. Also, requiring 
Cr{k) G iS(M) would be a simplification of the condition 2. However, this 
requirement for all Cr{k) might be too restrictive. 

Under the above assumptions we have 



where {g^xij))) is given by 



{9{x{t))) = g{x)h{x,r)dx. 



From assumption 2 it follows that Cr{k) G C 



N-li 



(3.5) 



(3.6) 



and thus we can 



conclude that for f{x) = e 



ikx 



^^ df{Z{T))dS{Z{T)) ^ ^ , df{x{T))dS{x{T)) ^ 

dZij) dZij) dxij) dx^r) 

J'fiZir))^ A-'fixir))^ 



Ei 



dZ^r) 



-) = {- 



dx'^^r) 



(3.7) 
(3.8) 
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and that the surface terms in the above integral expressions vanish. Thus 
h{x, t) obeys the pseudo F-P equation with a complex drift term 



dh{x, t) 
dr 



ld_ 

2dx 



dSix) d 
dx dx 



h{x, t) = Th{x, r). 



The operator Q, which is T transformed to L^(M, djl) by 

Q = e^f e"^, 

where djl is given by 

dfl{x) 



-S{x) 



Af 

has a zero eigenvalue with eigenfunction 

h{x) = 1. 



dx. 



(3.9) 

(3.10) 
(3.11) 

(3.12) 



However, we also note that Q has a second solution of zero eigenvalue in 
(M, djl) given by 

h jx) = r e^'^y^y . (3.13) 



Although e~^^^'^hs{x) is in general not positive definite, and hence auto- 
matically excluded when S is real, it is not a priori clear that /i(x, r) — >■ 
e~^^^%g[x) a.e. is excluded for complex 5*. 
But note 



,x 



0{- 



,N-1 ■ 



for |a;| oo 



(3.14) 



where N is the degree of the polynomial S{x). Thus e ^^^''hs^x) ^ S(M) , 
which contradicts |3.4| . So e~^^^^hs(x) cannot be the limit of h{x, r). As was 



true with the real S case, the finite measure /i implies that all solutions of 



3^ satisfying |3]^ converge pointwise to the desired result. 

The above 3 conditions appear to be a necessary set which must be satis- 
fied in order to connect the Langevin process Z to the complex Fokker Planck 
equation \i.9\ . These conditions also guarantee the correct convergence of the 
first A^ — 1 moments. More generally, however, any higher moment E{Z"'{r)) 
which does exist will converge correctly if h{x, r) is such that |(a;"(r))| < cxd. 

Two conclusions from this analysis are especially worth noting. First 
of all, we have seen that the complex Fokker-Planck equation has been de- 
rived under certain assumptions about the Langevin process. If these criteria 
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are not met, in particular if the function Cr{k) is ill behaved, then there is 
no longer necessarily a relation between the CL equation and the complex 
Fokker-Planck equation. This is true even if the spectrum of the complex 
F-P equation is such that all solutions converge to the desired stationary 
solution. Secondly, our result implies that the success or failure of the CL 
method only depends on the properties of the two eigenf unctions of the zero 
eigenvalue of Q and not on the requirement that the real part of the spec- 
trum of Q, Rea{Q) < . Thus, we see that, although it is certainly more 
convenient if Iiea{Q) < 0, an analysis of the spectrum is neither necessary 
nor sufficient for studying the behavior of the CL process. 



4 An Exactly Solvable Example 

For the Gaussian model S{x) = ux^ , u = a + ib E C^, the two dimensional 
real Fokker-Planck equation has the correct unique stationary density to 
which all initial solutions converge 0. The transition density is given by 

/(x,y,r|a;o,2/o,0)~e-^"W^-'W^M, (4.1) 

with 

ri{T) = X - mi{xo,yo,T) , r^^r) = y ~ m2{xQ,yQ,T) . (4.2) 

For the detailed form of the matrix S(r) and the vector m{xQ,yQ,T), the 
reader is referred to reference As mentioned above, the transition density 
converges to the stationary density 

Jim /(x,?/,r|xo,?/o,0) ~ exp[-aa;^ + 2^x?/ + (1 + ^)y^] . (4.3) 

With the transition density we have 

Crik) = E{e'''^^^^) = f e'^'dxdy [ f{x,y,T\xo,yo,0)f{xo,yo,0)dxodyo , 

(4.4) 

from which it follows that for /(xq, yo, 0) = 6{xo — x )S{yo — y) , Cr{k) has 
the form 

Cr{k) ~ e-'^'"^"'''") . (4.5) 

By examining the form of the matrix S(r), one can see that there always 
exists a tq such that for all r > tq, the real valued function a(a,6, r) > 0. 
Clearly, Cr{k) G 5(M), and the conditions 1 to 3 hold. 
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Since in this case the ground state of T is unique, by the preceeding 
analysis it follows that all initial states converge to the proper probability 
density. 

It is interesting to note that for quadratic polynomials it can be demon- 
strated that all solutions to the complex Fokker-Planck equation which are in 
L^(]R, dp,) converge to the ground state of Q or to zero. Note that all square 
integrable solutions are also absolutely integrable (L^) in this case. 

In this example the complex Fokker-Planck equation can be solved ex- 
actly. For the above range of the complex coupling u, the spectrum of Q 
satisfies Re[(j(Q)] < , and the eigenfunctions of Q defined on L^(]R, rf/t) 
form a basis. 

To see this note that 

Q = \e"\-H)e-''^ , (4.6) 

where H is now the ordinary Schrodinger operator for the complex frequency 
extension of the harmonic oscillator with a zero energy groundstate. With 
the help of the dilatation operator, which is given by ||T^ 

Ueij{x) = e''^i){e'x) , (4.7) 

one can now show the completeness of the eigenfunctions of H . First note the 
domain D{Ug) is dense in L^(R, dx). Secondly, the eigenfunctions of the real 
harmonic oscillator are in D{Ug) for Re[exp(2^)] > 0. The dilatation operator 
maps the eigenfunctions of the real harmonic oscillator to the eigenfunctions 
of the complex harmonic oscillator. Now for all functions 0(x) G D{Ug) and 
ipn{x) the n-th eigenfunction of the real harmonic oscillator we have 

(0, Ue^n) = (f/,^ V, ^n) = , for all n (4.8) 

if and only if U^^cj) = 0, which implies = 0. But since U^^cf) for G D{Ue) is 
dense, the eigenfunctions of the complex harmonic oscillator form a complete 
basis. Thus 

s — lim h(x,T) = hix) or (4.9) 
for all solutions h{x,T) G L^(M, d/i) of 

?^^ = Qh{x,r). (4.10) 
Again, from /i(M) < oo one can conclude that h{x, r) converges pointwise. 
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5 Conclusions 



To date there is no comprehensive theory of the complex Langevin method. 

General results are difficult to prove because many theorems on differential 
operators and diffusion processes do not apply to T and the Langevin equa- 
tions. However, for the case of polynomial actions we have demonstrated a 
set of sufficient conditions to guarantee convergence of the CL method to the 
correct results. 

An extension of these results from polynomial actions to the more general 
case of an analytic S requires a more detailed understanding of the solutions 
f{x,y,T) of T. Also of interest would be similar results for compact mani- 
folds. Work is currently progressing in these directions. 
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